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ABSTRACT 



Aims. We re-examine the effects of Compton scattering on the emergent spectra of hot DA white dwarfs in the soft X-ray range. Earlier 
studies have implied that sensitive X-ray observations at wavelengths A < 50 A might be capable of probing the flux deficits predicted by the 
redistribution of electron-scattered X-ray photons toward longer wavelengths. 

Methods. We adopt two independent numerical approaches to the inclusion of Compton scattering in the computation of pure hydrogen atmo- 
spheres in hydrostatic equilibrium. One employs the Kompaneets diffusion approximation formalism, while the other uses the cross-sections 
and redistribution functions of Guilbert. Models and emergent spectra are computed for stellar parameters representative of HZ 43 and Sirius B, 
. and for models with an effective temperature T c ff = 100000 K. 
(^) Results. The differences between emergent spectra computed for Compton and Thomson scattering cases are completely negligible in the case 
\Q of both HZ 43 and Sirius B models, and are also negligible for all practical purposes for models with temperatures as high as T eff = 100 000 K. 
Models of the soft X-ray flux from these stars are instead dominated by uncertainties in their fundamental parameters. 
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1. Introduction In a later paper. iMadel dl998b computed the effects of 

Compton scattering for a model corresponding to the pa- 

. £ ; The scattering of radiation by free electrons is one of the dom- ramet ers of the DA white dwarf HZ 43. Differences between 

, inant sources of continuous opacity in the atmospheres of hot Compton and Thomson scattering model spectra were apparent 

in white dwarfs. Most model atmosphere calculations adopt the for ^ < 50 A, and grew to orders of magnitude by 40 A. While 

.«? classical Thomson isotropic scattering approach, whereby only current x . ray instrume ntation is not sufficiently sensitive to 

the direction of photon propagation changes as the result of a study the spectra of even the brightest D A white dwarfs in any 

scattering event. More rigorously, finite electron mass instead detail at wave l e ngths A < 50 A, the spectral differences im- 

implies that both momentum and energy exchange should ac- plied by the more rigorous Compton redistribution formalism 

tually occur. w jjj be Q f m t ereS ]; t m0 re sensitive future missions. Moreover, 

The effects of Compton scattering on whit e dwar f mode l the Chandra Low Energy Transmission Grating Spectrometer 

atmospheres was first investigated in detail bv lMadeil dl994 . (LETG+HRC-S) effective area calibration is based on observed 

who found that pure hydrogen models with temperatures of spectra of HZ 43 and Sirius B at wavelengths A > 60 A 

10 5 K show a significant depression of the X-ray continuum JPease etal1l2003h . It is therefore of current topical interest 

for wavelengths < 50 A. Effects for models containing signifi- to re-examine the influence of Compton scattering for these 

cant amounts of helium, or helium and heavier elements, were stars and determine whether any significant differences might 

found to be much smaller or negligible, in keeping with expec- be discernible between Thomson and Compton scattering in the 

tations based on the relative importance of electron scattering LETGS bandpass, 
as an opacity source as opposed to photoelectric absorption. 

In this paper, we perform two independent and rigorous 

Send offprint requests to: V. Suleimanov tests of the influence of Compton Scattering on the emergent 

Correspondence to: e-mail: suleimanov@astro.uni-tuebingen.de spectra of hot DA white dwarfs. Our methods of calculation 
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are outlined in |2 while results and conclusions are briefly dis- 
cussed in §^3]- 13 

2. Computational Methods 

In both our numerical approaches, outlined below, we com- 
puted model atmospheres of hot white dwarfs subject to the 
constraints of hydrostatic and radiative equilib rium assuming 
planar geometry using standard methods (e.g. Mihalas 1978). 
The equation of state of an ideal gas used assumes local ther- 
modynamic equilibrium (LTE), and therefore did not include 
terms describing the local radiation field. 

The model atmosphere structure for a hot WD is described 
by the hydrostatic equilibrium equation, 



dP„ GMy, 



dm R 2 , 

wd 



- 47T 



Hy dv, 



(1) 



where k v is opacity per unit mass due to free-free, bound- 
free and bound-bound transitions, <x v is the electron (Compton) 
opacity, H v is Eddington flux, P g is a gas pressure, and m is col- 
umn density 



dm 



-pdz . 



(2) 



Variable p denotes the gas density and z is the vertical distance. 
As is obvious from Eqn. ^ the structure of the atmosphere is 
coupled to the radiation field and the structure and radiative 
transfer equations need to be solved simultaneously under the 
constraint of radiative equilibrium. 

In the Thomson approximation, in which no energy or mo- 
mentum between photons and electrons is exchanged, cr v — cr e , 
where <x e is the classical Thomson opacity. 

2.1. Method 1 

In our first approach, Compton scattering is taken into account 
in the radiat ion transf e r equa ti on using the Kom p aneets 
operator jKompaneetsI 1 19571 IZavlin & Shibanovl Il99ll 
Greb enev & Sunvaevl2002l) : 



d fvJv _ k y 
dr}. k v + (r. 



(J V ~ By) 



cr e kT 
ky + cr e m»c 2 



(3) 



d ( dJy Teff / CJy\\ 



where x = hv/kT e g is the dimensionless frequency, / v (r v ) « 
1 /3 is the variable Eddington factor, J v is the mean intensity of 
radiation, B v is the black body (Planck) intensity, T is the local 
electron temperature, T e s is the effective temperature of WD, 
and C = c 2 h 2 / 2(kT e ff) 3 . The optical depth t y is defined as 



dT v = {k v + cr e ) dm. 



(4) 



These equations have to be completed by the energy balance 
equation 



f 

Jo 



ky (Jy ~ By) flfV ~ (T S 



kT 

m e c 2 



(5) 



^ jT JydV-^f f~ Xjy(l + rfvj = 0, 



the ideal gas law 

P s = N ln ,kT, 



(6) 



where iV tot is the number density of all particles, and also by 
the particle and charge conservation equations. We assume lo- 
cal thermodynamical equilibrium (LTE) in our calculations, so 
the number densities of all ionisation and excitation states of 
all elements have been calculated using Boltzmann and Saha 
equations. 

For solving the above equations and computing the model 
atmosphere we used a version of the computer code ATLAS 
dKuruczll 970Hl993 h. mod ified to deal with high tem peratures; 
see llbragimov et alTd2003l) and Swart z^t alJ J2002) for further 
details. This code was also modified to account for Compton 
scattering. 

The scheme of calculations is as follows. First of all, the 
input parameters of the WD are defined: the effective temper- 
ature r e ff and surface gravity g = GM wd //?^ d . Then a starting 
model using a grey temperature distribution is calculated. The 
calculations are performed with a set of 98 depth points m\ dis- 
tributed logarithmically in equal steps from m ss 1CT 6 g cirT 2 to 
'«max- The appropriate value of m max is found from the condi- 
tion -\/r Vi b-f,f-f(?w max )T v (m max ) > 1 at all frequencies. Satisfying 
this equation is necessary for the inner boundary condition of 
the radiation transfer. 

For the starting model, all number densities and opac- 
ities at all depth points and all frequencies (we use 
300 logarithmically equidistant frequency points) are cal- 
culated. The radiation transfer equation (0 is non-linear 
and i s solved iteratively by the Feautri e r method dMihalasI 



19781 see also IZavlin & Shibanovl [l99ll IPavlov et all 1 199 it 
Gre benev & Sunva ev 2002). We use the last term of the equa- 
tion (0 in the form xJ' v (l + CJ' V X /x 3 ), where J' v ~ l is the mean 
intensity from the previous iteration. During the first iteration 
we take /Jr 1 = 0. Between iterations we calculate the variable 
Eddington factors f v and h v , using the formal solution of the 
radiation transfer equation in three angles at each frequency. 
Usually 2-3 iterations are sufficient to achieve convergence. 
We used the usual condition at the outer boundary 

where h v « 1 /2 is surface variable Eddington factor. The inner 
boundary condition is 

dJy _ dB 



(7) 



(8) 



The outer boundary condition is found from the lack of incom- 
ing radiation at the WD surface, and the inner boundary condi- 
tion is obtained from the diffusion approximation J v « B v and 

Hy * \/3XdBy/dTy. 

The boundary conditions along the frequency axis are 



Jy — By 

at the lower frequency boundary (v ff 
kT e ff) and 



dJy - . Teff / CJy 

X— - 3Jy + —Xjy II + — 







(9) 

10 12 Hz, hVmin <K 



(10) 
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at the higher frequency boundary (v max * 3 • 10 17 Hz, hv max » 
kT e ff). Condition (|9jl means that at the lowest energies the true 
opacity dominates the scattering k v » <x e , and therefore J v ss 
B Y . Condition dlOt means that there is no photon flux along 
frequency axis at the highest energy. 

The solution of the radiative transfer equation Q was 
checked for the energy balance equation together with the 
surface flux condition 



An 



H v (m = Q)dv = oT e 4 ff = 4nH 



The relative flux error as a function of depth, 
s H (m) = 1 - 



H v (m)dv 



(11) 



(12) 



was calculated, where H v (m) is radiation flux at given depth. 
This latter quantity was found from the first moment of the 
radiation transfer equation: 

dfyJy 



djy 



= Hy 



(13) 



Temperature corrections were then evaluated using three dif- 
ferent procedures. The first is the integral A-iteration method, 
modified for Compton scattering, based on the energy balance 
equation Q. It is valid in the upper atmospheric layers. The 
second one is the Avrett-Rrook flux correction, which uses the 
relative flux error, and is valid in the deep layers. And the third 
one is the sur face corr ection, which is based on the emergent 
flux error. See iKuruczl \ 1 970|) for a detailed description of the 
methods. 

The iteration procedure is repeated until the relative flux er- 
ror is smaller than 1%, and the relative flux derivative error is 
smaller than 0.01%. As a result of these calculations, we ob- 
tain the self-consistent WD model atmosphere together with 
the emergent spectrum of radiation. 

Our method of calculation was tested on a model of burst - 

)l HMadeill99ll) . 



ing neutron star atmospheres ( Pavlo v et alJH 
and a model DA wh ite dwarf atmosphere with — 2 • 10 5 K, 
log g = 6.3 llMadeill994 . A greement with the earlier calcula- 
tions is extremely good. We show the emergent spectrum from 
the latter calculation in Fig^ 

2.2. Method 2 

Our second approach adopts the equat ion of transfer fo r ab- 
sorption and scat tering presented by Sampson] Jl959l) and 
IPomraningl dl973l see Eqn. (2.167)). The equation of transfer 
can be expressed in the form 

d^fyJy ky 



k v + cr. 



(Jv ~ By) + 



oo 

^—Jy \ <£(V,V0 (l 

+ cr v J \ 



2hv' 



Jy> dV 



ky + (Ty 



1 + 



2hv 3 



Jy \ X 



j ^(v.vVv'^j exp 



h(v - Y) 



kT 



dv . 



(14) 
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Fig. 1. Spectra (a) and temperature structures (b) of a hot 
DA white dwarf model atmosphere with r e jf = 200000 K 
and log g = 6.3 with (solid line) and without (dashed line) 
Compton scattering taken into account. 



This transfer equation is written on the monochromatic opti- 
cal depth scale dr v = —(k v + <Ty)pdz. The variable J v denotes 
the energy-dependent mean intensity of radiation. The func- 
tion <t>(v, v') denotes the zeroth angular moment of the angle- 
dependent differential cross-section, normalized to unity (see 
below). 

Tr ansformation of the equation of transfer by IPomraning l 
Jl973li to the Eq. Eland the definition of required angular ap- 
proximations f or Compton scatterin g i n a stellar atmosphe re 
was outlined bv lMadeil ( 1 1 99 ill 1 994 and lMadei et alJ J2004 . 

Our equations and theoretical models of Method 2 use 
detailed differential cross-sections for Co mpton scattering , 
cr(v — > v',n ■ n'), which were taken from iGuilbern 41981). 
Cross-sections correspond to scattering in a gas of free elec- 
trons with relativistic thermal velocities, and they are also com- 
pletely valid at low temperatures. Differential cross-sections 
were then integrated numerically to obtain large grids of 
Compton scattering opacity coefficients cr v 



dv' cr(v — » Y , n ■ n') , 



(15) 
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and grids of angle-averaged Compton scattering redistribution 
functions (D(y, v') 



1 C doj f 
0(v, V) = — Q) — cr(v -> v , n ■ ri) . 



(16) 



The scattering frequency redistribution function was intro- 
duced by Po mraniiigl(ll973l) — see Eqn. (7.95) of his book — and 
it represents the normalized probability density of scattering 
from a given frequency v to the outgoing frequency V . 

Note that Eq. El also includes stimulated scattering terms, 
1 + (c 2 /2hv y )J v , which ensure the physically correct descrip- 
tion of Compton scattering. The actual equatiqns_and calcu- 
lations used here strongly differ from those in Madei ( 1998). 
The physics of Compton scattering used here is also funda- 
mentally different to Method 1 . Note, that the latter method and 
that of Madei ( 1998) employ the well-known Kompaneets dif- 
fusion approximation to Compton scattering kernels. The ap- 
parent difference between Methods 1 and 2 is that they use ei- 
ther the differential Kompaneets kernel (Method 1) or kernels 
given by integrals over the detailed Compton scattering profiles 
(Method 2). 

The equation of radiative equilibrium (the energy balance 
equation) requires that 



f 

Jo 



Hydv 



4n 



(17) 



The above condition is fulfilled in a hot stellar atmosphere, 
where energy transport by convective motions can be ne- 
glected. 

Computing derivatives of both sides of Eqn. ^] and using 
the equation of transfer, Eqn. [2] one can obtain the alternative 
energy balance equation 



OO OO OO 

j ky(Jy-By)dV+ ^(TyJydV ^0^,/)^ 



CO 

-M 



2hv' 3 



Jy> dV 



J J y ®{v,V){^J exp 



h(v - /) 



kT 



dV=0, (18) 



which can be compared with its analog, the Eqn.[5]of Method 
1. 

The computer code atm2 1 used for the model cal cula- 
tions was described in detail in Madei & Rdzariska (2000) and 
Mad ei etaD J2004I) . The struct ure of the code is based on the 
partial linearization scheme by Mihalas ( 1978), in which cor- 
rections of temperature AT" and the function AO are built into 
the equation of transfer. The high numerical accuracy and very 
good convergence properties of the atm21 code are vital for 
the present paper, and allowed us to compute accurate spectra 
for atmospheric parameters appropriate for the white dwarfs 
HZ 43 and Sirius B using Method 2, outlined above. These cal- 
culations supersede less accurate X -ray spectra o f HZ 43 which 
were presented in the earlier paper (Mad eill998l) . 



For the present research, the atm21 code included numer- 
ous bound-free LTE opacities of neutral hydrogen and free-free 
opacity of ionized hydrogen, which always remains in LTE. 
No hydrogen lines were included in the actual computations. 
In each temperature iteration the code solves the equation of 
hydrostatic equilibrium to obtain stratifications of gas pressure 
P g and density p in the model atmosphere. After that the atm21 
code solves the set of coupled equations of radiative transfer 
with implicit temperature corrections and finds the stratifica- 
tion of AT in the model atmosphere. The equation of radiative 
transfer was solved using the Feautrier method and the tech- 
nique of variable Eddington factors (Mih alasl 19781) . Boundary 
conditions along the t v axis were the same as in Method 1. 
Explicit expressi ons for the tempera ture corrections, AT, can 
be found, e.g., in iMadei et alJJ2004 . 

3. Computations and Results 

Using our two independent methods, we investigated the ef- 
fects of Compton scattering on the emergent spectra of three 
pure H white dwarf models. These were: models appropri- 
ate for the well-known DA white dwarfs Sirius B (assuming 
r eff = 24700 K, log g = 8.6) and HZ 43 (r eff = 50000 K, 
log g = 8.0), together with a significantly hotter model with 
T e ff = 1 • 10 5 K and log g = 6 and 8. These models cover the 
effective temperature range relevant to pure H DA atmospheres 
and specifically address the question of whether Compton scat- 
tering might be relevant to the X-ray bright DA white dwarfs 
Sirius B and HZ 43 that are central to the low-energy calibra- 
tion of Chandra JPease et all20 03). 

The results of our calculations from both methods are pre- 
sented in Figures l2"l5l In Fig.|3^ the spectra of the model atmo- 
sphere for the DA white dwarf HZ 43 (r eff = 50 000 K, log g = 
8.0) computed using Method 1 with (solid line) and without 
(dashed line) Compton effects are shown. We also calculated 
a non-LTE model atmosphere for H Z43 using the Tubingen 
Model Atmosphere Package (TMAP) JWerner et all2003l) and 
computed the radiation transfer equation Q with Compton 
scattering using this non-LTE model atmosphere structure. The 
corresponding spectra are shown in Fig.|2k by the dotted line 
and by open circles. 

Calculations using Method 2 in LTE are illustrated in Fig. [3] 
The differences between Compton and Thomson models in 
these calculations are very similar to the differences found in 
Method 1 in Fig.|2t, and are significant only at wavelengths < 
25 A. The flux in the models with Compton scattering (Method 
1) is smaller than the flux in the model without Compton scat- 
tering by 0.5% at 50 A, 1.3% at 40 A, 7% at 30 A, and 40% at 
20 A. Similarly, model computations performed with Method 2 
yield the following differences for HZ 43: 0.5% at 50 A, 1.4 % 
at 40 A, 5% at 30 A, and 56% at 20 A. 

These results can be understood more clearly if we consider 
the Comptonisation parameter Zc om pt- 

hv 



Zc 



ompt 



■max((r*) ,t*), 



(19) 



where hv/m e c 2 is the relative photon energy lost during 
one scattering event off a cool electron, max((r*) 2 , t*) is 
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Fig. 2. Spectra (a) of the DA white dwarf HZ 43 model at- 
mosphere with (solid line) and without (dashed line) Compton 
scattering using Method 1. Also shown are the spectra of the 
non-LTE model atmosphere (dotted line) and the spectra of 
non-LTE model atmosphere with Compton scattering (open 
circles). Run (b) of Comptonisation parameter Zc om pt with 
wavelength (see Eqn.lT^. 



the number of scattering events the photon undergoes be- 
fore escaping, t* is the Thomson optical depth, correspond- 
ing to the depth where escaping photons of a given fre- 
quency are created. The Comptonisation parameter is, then, 
a representation of the influence of Compton down-scattering 
on the emergent spectrum: significant Compton effects are 
expected if the Comptonisa tion parameter approaches unity 
jRvbicki & Lig htma nT ll979l) . In FiglJJ) the dependence of 
Zcompt on the wavelength is shown. It is clear that the 
Comptonisation parameter is very small down to 25 A, and this 
is indeed reflected in the emergent spectrum. 

Similar results were obtained for the DA white dwarf 
Sirius B (assuming T e ff = 24700 K, log g = 8.6). The corre- 
sponding spectra and the Comptonisation parameter are shown 
in Fig. [4] The effective temperature of this star is smaller, and 
the surface gravity is larger, therefore, Compton scattering is 
even less significant than for HZ 43. We also calculated a non- 
LTE model atmosphere for Sirius B using the TMAP and com- 
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Fig. 3. Emergent X-ray spectra for a pure hydrogen model at- 
mosphere appropriate for the parameters of HZ 43 computed 
for classical Thomson scattering (dashed) and C ompton scat- 
tering (solid) using Method 2 (the method of iMadei et alJ 
2004). 



puted the radiation transfer equation (0 with Compton scat- 
tering using this non-LTE model atmosphere structure, as for 
HZ 43. The corresponding spectra are shown in Fig.|4k by the 
dotted line and by open circles. Compton scattering changes 
only very slightly the emergent spectrum at wavelengths below 
20 A. 

In the Fig. |3 we present the spectra of the hot DA white 
dwarfs model atmospheres (T e tf = 100000 K, log g — 6 
and 8) with and without Compton scattering. It is obvious that 
Compton effects are more significant for these hot DA mod- 
els (especially with low surface gravity) than for Sirius B and 
HZ 43, but visible effects still do not occur for wavelengths 
> 50 A. Of course, log g — 6 is not a so realistic value for 
white dwarfs, but we calculate this model to demonstrate the 
dependence of the Compton effect on the surface gravity. 

4. Discussion 

The differences we find here between emergent spectra com- 
puted using Thomson and Compton scattering are essentially 
negligible for practical purposes. Such differences are in fact 
much smaller than the predicted spectral uncertainties result- 
ing from uncertainties in the current knowledge of the funda- 
mental parameters of DA white dwarfs - even the best known 
examples such as HZ 43 and Sirius B. The spectrum in this 
Wien tail region is especially sensitive to uncertainties in the 
effective temperature. 

Comparison of the calculations present ed here using both 
Methods 1 and 2 with the earlier work of lMadeil dl998l) do, 
how ever, reveal s ome significant differences. For HZ 43, Fig. 3 
of iMadeil ( 119981) suggests a large X-ray flux deficit due to 
Compton scattering of photons to longer wavelengths, with a 
precipitous decline in emergent flux at ~ 40 A. However, our 
Figs.|2]and[3]illustrate much small er effects. 

A careful examination of the IMadeil dl998l) computer code 
performed by one of us (JM) has shown that the earlier code 
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Fig. 4. The same as Fig. [2] but here for the DA white dwarf 
Sirius B LTE model atmosphere. The dotted line lies every- 
where under the open circules (except for hydrogen lines). 



(which used Kompaneets scattering terms) could strongly ex- 
aggerate effects of Compton scattering in cases when they were 
of only marginal significance. This is the case for the X-ray 
spectrum of HZ 43. The effect was caused by an approxima- 
tion adopted in the solution of the transfer equation that was 
quite valid for the study of X-ray burst sources, for which the 
code was primarily developed, but which became marginally 
inaccurate for the case of hot DA white dwarf atmospheres. 
This problem has been solved by the very stable algorith m of 
the new code (see Method 2), described in Madei & Rozanska 
i2000fc . 

One sho uld note that t he differences between this work 
and that of Madei (1998) are not related to the numeri- 
cal approaches adopted for Compton scattering. Both the 
Kompaneets diffusion approximation (Method 1) and the 
Compton scattering terms of the integral form (Method 2) satis- 
factorily describe effects of Compton scattering in X-ray spec- 
tra of hot white dwarf stars. The results of the work prese nted 
here using both methods supersede those of Madei Jl998l) . 

At the higher effective temperatures represented by the 
r e ff = 100000 K models, the emergent spectra for Thomson 
and Compton scattering begin to diverge at ~ 50 A for the 



Fig. 5. Spectra (a) of the DA white dwarf r eff = 100000 
K model atmospheres with and without Compton scattering. 
Solid lines correspond to model with log g = 8.0, and dashed 
lines correspond to model with log g = 6.0. Softer high en- 
ergy tails correspond to models with Compton effect. Run (b) 
of Comptonisation parameter Zc om pt with wavelength. 



log g — 6 model, and the effects are much larger toward shorter 
wavelengths than for the higher gravity logg = 8 case. It is 
of course questionable as to whether any white dwarfs with 
pure H atmospheres exist with such high effective tempera- 
tures, since in hotter stars radiative levitation tends to enrich the 
atmosphere with metals. In atmospheres with significant metal 
abundances the electron scattering opacity is insignificant com- 
pared with that due to metals, and Compton redistribution ef- 
fects are rendered irrelevant. 

Our calculations show that even in the hottest pure H at- 
mospheres it is highly unlikely that future X-ray observations 
will be sufficiently sensitive to discern the Compton redistri- 
bution effects. Again, uncertainties in model parameters such 
as effective temperature and surface gravity, and abundances of 
He and trace elements, together with uncertainties in param- 
eters entering into the modeling calculations themselves will 
dominate. Even for models normalised to the same flux at UV 
wavelengths a 1% error in r e tf will induces a 20% flux error at 
75A. 
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5. Conclusions 

New calculations using two independent, rigorous numerical 
methods confirm that the effects of Compton energy redis- 
tribution in photon-electron scattering events are completely 
negligible for the interpretation of X-ray spectra of DA white 
dwarfs such as Sirius B and HZ 43. Differences between emer- 
gent spectra of Compton and Thomson cases are in fact much 
smaller than the predicted spectral uncertainties resulting from 
uncertainties in the current knowledge of the fundamental pa- 
rameters of stars such as HZ 43. We have found that differ- 
ences between Compton effects predicted here for HZ 43 and 
the calculations of Madei ( 1998) are caused by approximations 
used in the solution of the transfer equation in the former work; 
the results presented here supersede the earlier ones. We con- 
clude that current non-LTE model atmosphere spectra of hot 
DA white dwarfs neglecting Compton scattering can be safely 
used for the calibration low energy detectors of X-ray observa- 
tories for the foreseeable future. 
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